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I. INTRODUCTION 

In recent paper [l] Rakos and Schiitz have solved the master equation for the one- 
dimensional totally asymmetric exclusion process (TASEP) with particle-dependent rates 
defined on the integer lattice Z. The TASEP with particle-dependent rates |2| generalizes 
the usual TASEP, where particles hop independently in continuous time with equal rates to 
their right neighboring sites provided the target sites are empty DO . A random choice of 
rates leads to condensation phenomena which have been studied intensively for the associated 
Z ero rang e processes QflfiU fl Q Q. Specific^, the m c,ion of the lowest p artic ,e 
generates a macroscopic traffic jam viewed in the coarse-grained scale as a discontinuity in 
the particle density. The solution of the master equation for the particle-dependent rates in 
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jj| was obtained under the general approach developed by Schiitz 0, 0] for the usual TASEP 
for a P-particle system on the infinite one-dimensional lattice. The conditional probability 
to find P particles in positions X\, X2, ■ ■ ■ , ip; x\ < X2 < ■ ■ ■ < xp at continuous time t if 
these are in positions x®, x%, . . . , x P ; x\ < x\ < ... < x° P at initial moment of time was found 
in the form of a determinant of P x P matrix. In this paper, we extend the determinant 
formula for more complicated case, the TASEP in discrete time with backward sequential 
update on the lattice with periodic boundary conditions. The usual TASEP with the ring 
geometry has been considered in 

0, Q- 

Two methods, an analytical Bethe Ansatz 
and a combinatorial solution were presented. Here, we use the combinatorial approach jlfil ]. 
based on a geometrical interpretation of the Bethe ansatz. 

The discrete time TASEP with backward-sequential update is defined as follows. At 
each discrete moment of time T, particles are tested for mobility successively in backward 
direction starting from an arbitrary empty site. A particle is movable at the moment T 
if its right site is not occupied by next particle at time T + 1. (In the parallel update all 
particles are tested simultaneously and the particle is movable at time T if its right site is 
empty at time T). Using a discrete time formulation, we interpret the Bethe ansatz as a 
combinatorial problem of cancellation of forbidden space-time trajectories of particles. A 
resulting expression for the conditional probability is then transformed into a continuous- 
time formula by changing the one-particle discrete Bernoulli dynamics for the continuous 
Poisson process. 

The paper organized as follows. The content of section II is a description of the com- 
binatorial approach, which is included for convenience and completeness. In subsection 
A of section II, we describe a combinatorial ansatz for directed walks using as an exam- 
ple the free-fermion problem of "vicious" walkers In subsection B, we introduce the 
combinatorial treatment of the Bethe ansatz and consider the discrete-time TASEP with 
same hopping probability for all particles. The main subject of the paper, the TASEP with 
particle-dependent rates is considered in section III. In section IV, we close the paper with 
some conclusions. 
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II. THE COMBINATORIAL METHOD 
A. Combinatorial ansatz for directed walks 

Consider P particles labelled {1, 2, ... , P} hopping in one direction on an infinite one- 
dimensional lattice. The discrete space-time dynamics of their motion can be described 
as a set of continuous broken trajectories on a triangle lattice A, which is obtained from 
the square lattice by adding a diagonal bond between the upper left corner and the lower 
right corner of each elementary square. Let (x, T) be integer space-time coordinates of 
a particle on A, where the vertical time axis is directed down and the horizontal space 
axis is directed right. A trajectory of each particle is a sequence of connected vertical and 
diagonal bonds of A. The diagonal bonds correspond to jumps of particles to their right 
for a unit time-step with probability v. The vertical bonds correspond to stays of particles 
at fixed sites during the unit time-step with probability 1 — v. The initial positions of 
particles are x\ < x® < • • • < x° P . The problem of "vicious" walkers is to find the 
conditional probability Pt{x\xP) for the particles to reach the positions x = (x\, X2, ■ ■ ■ , xp); 
x\ < X2 < ■ ■ ■ < xp from the initial positions x° = (x^x®, . . . ,x° P ) for time T so that no 
pairs of trajectories are intersected during time T. It means that every site can be occupied 
by at most one particle at every moment of discrete time and if this rule is violated at a 
moment T' < T, the process stops. 

We start with consideration of one-particle motion on the infinite lattice. Let Tt(x\x°) be 
a set of one-particle trajectories, which are starting at (x°, 0) and finishing at (x, T). Each 
trajectory q e Tt{x\x°) is realized with probability v x ~ x °(l — v) T ~ x+x ° . Therefore, the total 
probability for the particle to reach x from x° for time T is 

P T (x\x°) = v x ~ x \l - *;) t -* +:e0 ||T t (x|x )|| = F (x - x°\T). (1) 

where 

F (x\T) = (^jv x (l-v) T - x (2) 

is the statistical weight of the set Tt(x\x°) of one-particle trajectories. For the case of P 
particles, the set Sp of all possible free trajectories is a direct product of one-particle sets of 
trajectories reaching (x,T) from (x°,0): 

§ P = 7^(xi|x?) <8> T T (x 2 \x° 2 ) <g> . . . <g> T T (x P \x P ) (3) 
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The set Sp contains non-intersecting and intersecting trajectories. The latter cases are un- 
allowed by the single-occupation rule and should be subtracted in evaluations of probability 
P T (f|f°). 

Consider first the case P = 2. We denote the set of unallowed elements by U2 C §2- 
To cancel the contribution of forbidden elements, we introduce an auxiliary set of pairs of 
trajectories 

A 2 i = T T (x 2 \x° l )®T T (x 1 \x° 2 ), (4) 

where the final coordinates are permuted with respect to those in § 2 - It is easy to see that 
all elements of the set A21 are pairs of crossing trajectories. Each intersecting pair q e 82 has 
a first collision point (x c ,T c ), i.e. the space-time point, where the trajectories meet for the 
first time. Consider the set §2(^0 T c ) of intersecting trajectories with the fixed first collision 
point (x c ,T c ). The sets § 2 (a; c , T c ), (x c ,T c ) 6 A break the set U 2 into subsets parameterized 
by coordinates (x c ,T c ). For all (a^,T c ') 7^ (x",T") we have 

§2«,T c ')n§ 2 «,T c ") = (5) 

and 

|J § 2 (x c ,T c ) = U 2 (6) 

(x c ,T c )eA 

For every set S 2 (x c , T c ), there exists a uniquely defined set A 2 i(x c ,T c ) C A 2 i obtained 
from §2(^0, T c ) by permutation of tails of all trajectories beginning at the first collision point. 
The sets §2(2^, T c ) and A 2 i(x c ,T c ) are geometrically equivalent, because there is one-to-one 
correspondence between their elements. Therefore ||§2(^ c , T c )\\ = ||A 2 i(x c , T c )\\. Like the 
sets §> 2 (x c ,T c ), the sets A 2 i(a; c ,T c ) break the set A 2 i into subsets: for all (x' c , T c ') 7^ (x",T") 
we have 

A 21 (x' c ,T c ')nA 21 «,T c ") = (7) 

and 

|J A 21 (x c ,T c ) = A 21 . (8) 

(x c ,T c )GA 

It follows from Eqs. (j5l6l7l8|) that the whole sets U2 and A 2 i are geometrically equivalent 
and HU2II = ||A 2 i||. Given (x,T) and (x°,0), all elements of §2 have the same weight 

2 

Q = Y[v x >- x °>{l-v) T - x * +x l (9) 

i=l 
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Then, we have for the probability 

P T (x\x°) = Q x (USall - ||A 21 ||) = (10) 
= (F (xi - x° 1 \T)F (x 2 - x° 2 \T) - F (x 2 - x° 1 \T)F {x 1 - x° 2 \T)) = 
= detM, 

where the elements of the 2x2 matrix M are 

M iJ = F (x j -x° i \T), i,j = 1,2. (11) 

This result can be easily generalized for the P-particle system of vicious walkers Q]. The 
set Sp contains the subset of intersecting trajectories Up. We introduce P\ — 1 auxiliary 

sets 

K(i2...p) = T T (x n{1) \xl) (8) T T (x n{2 ) \xl) ® . . . <8> T T (x n(P) \x° P ) (12) 

where 7r is any permutation of numbers 1,2, ... ,P beside the identical one 1. Each element 
of the set 

Up|J (\J \ K(i,2,..,P) J (13) 

containing trajectories intersecting at given point (x c ,T c ), has a unique geometrically 
identical counterpart with a pair of permuted indices, e.g. 1,2, . . . ,i, . . . , j, . . . , P and 
1,2, . . . , j, . . . ,i, . . . , P. Introducing the sign of permutation sgn(ir), we can ascribe opposite 
signs to two counterparts. Then, we obtain 

P T (f|x°) = Qx Lsp\\ +Y,s9n(*)\\K(i2...p) II j = (14) 



where 



and 



det M, 



M i;j = F ( Xj - x^T), i, 3 = 1,2,... P (15) 



p 

Q = Y[v Xi - x °(l-v) T ~ Xi+x l (16) 
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B. Combinatorial treatment of Bethe ansatz 

The vicious walkers are locally interacting particles. In this subsection, we consider the 
system of particles with the interaction of a non-zero range, which takes place in the case 
of the totally asymmetric exclusion process with the backward sequential update |3, Q . 
Consider again P particles hopping in one direction on an infinite one-dimensional lattice. 
The interaction between particles can be defined by the following rules: 

1. Trajectories of particles do not intersect (every site can be occupied by at most one 
particle. 

2. A particle stays at its own site with probability 1 if the target site is occupied by 
another particle during the step of discrete time. 

Like the previous case, the trajectory of each particle is a sequence of vertical and diagonal 
bonds on the lattice A. Each diagonal bond has weight v and the vertical bond weight 
1 — v. In view of more complicated interaction, it is convenient to decompose the set of 
all free trajectories into more detailed subsets. For each vertical bond [(x,t), (x,t + 1)], the 
trajectory of z-th particle passing this bond can be decomposed into two parts 0) — > 
(x, t)\l\(x, t + 1) — ► (xi,T)] and [(x°,0) — > (x,t)\ — v\(x,t+l) — ► (xi,T)] where the value 
between vertical bars means the weight of the bond [(x, t), (x,t + 1)]. These new trajectories 
are geometrically equivalent, but the first trajectory passes the selected vertical bond with 
weight 1 and second one with weight —v. We make this decomposition for each vertical bond 
of each trajectory and denote the whole set of decomposed trajectories of z-th particle by 
Tt(xi\x®) using the same notation as for the set of one-particle trajectories in the previous 
subsection. The weight of the trajectory q E T T ( product of the weights of its 

bonds. 

T 

fj,(q) = Y\ (i(k-th bond of q) (17) 

fc=i 

The weight of the set of trajectories is a sum of weights of its elements. In accordance with 
the above definition (JTJ), the weight of the whole set 7r(^j|x°) is 

fi(T T (x l \x° i )) = F (x t -x° i \T) (18) 

The set of free P-particle trajectories Sp is defined by Eq. Q and has the weight 

P 

MSp)= EM9) = n^(l-«) T ^(_ ^ n). (19) 

geSp i=i \ * tJ 
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As in the case of vicious walkers, the set Sp contains a subset of unallowed elements,Up, 
which should be excluded. By the rules 1 and 2, an element of S P is unallowed, if there is 
at least one pair of intersecting trajectories, or if there are two neighboring vertical bonds 
from which the left one has weight —v. To cancel unallowed elements we start as above with 
the case P = 2 and construct an auxiliary set of trajectories A21. We will show that, due to 
non-locality, the auxiliary set has more complicated structure 



fc 1= 0fc 2 =0 v 7 

Beginning the construction, we notice that every unallowed pair has a first collision point 
(x c , T c ) G A, where the particles come for the first time to neighboring sites at a moment 
T c < T. The first trajectory reaches site (x c , T c ) from (x?, 0) and the second one reaches site 
(x c + 1, T c ) from (x®, 0) and then it makes the vertical step to site (x c + 1, T c + 1). 

The first trajectory has, just after the collision, a diagonal bond [(x c , T c ), (x c + 1, T c + 1)] 
with weight v (referred to as a collision of the first type) or vertical bond [(x c , T c ), (x c , T c + 1)] 
with weight — v (refereed to as the collision of the second type). Notice, that trajectories, 
which have the vertical bond [(x c , T c ), (x C) T c + 1)] with weight 1 are allowed. Thus, we have 
two types of unallowed elements of U 2 for fixed (x c ,T c ). Denote these subsets by Y(x c ,T c ) 
for the first type and W(x c ,T c ) for the second one. For all (x' c ,T c ') 7^ (x",T") we have 



For each space-time point (x c , T c ) e A, we construct the sequences of pairs of trajectories 
A v (ki,k 2 ) and A w (ki,k 2 ) with k\ = 0,1,2,... and k 2 = 0,1 (arguments x c and T c are 
omitted), obtained from V(x c ,T c ) and W(x c , T c ) by permutation of tails in each pair of 
trajectories and shifting the initial parts in negative direction. Explicitly, the transformation 




(20) 



Y(x' c X) n V(x' c ',T c ")=0, 
W«,T C ') n W«,T c ")=0, 

WcX) n w«,t c ")=0, 



and 




(21) 
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V(x c ,T c ) A v (k u k 2 ) is 



( «0) \ 

I 

(x c , T c ) 
i 

(x c + l,T c + 1) 
I 



<8> 



( 04 o) ^\ 

(x c + l,T c ) 
I 

x c + l,T c + l) 
i 

V (^,t) y 



( (x?-ifei,0) \ 
I 

(x c fci, T c ) 
I 

x c + l-fci,T c +l) 
I 

(*2,r) y 



/ (x° - h, 0) 
I 

(x c + 1 - fc 2 , T c ) 
i 

x c + l-k 2 ,T c + l) 
i 

\ (xi,T) I 



and the transformation W(a; c ,T c ) =>- K w (k\, k 2 ) is 



/ (*?,0) \ 
I 

(x c , T^) 
I 

(x c ,T c +l) 
I 



( 04 0) \ 
I 

(x c + 1,T C ) 
I 

x c + l,T c + l) 
I 

V (^,t) y 



( (xS-fci,0) \ 

I 

I 

(x c -fci,T c + l) 
I 

V fe,T) y 



(x° 2 - k 2 , 0) \ 
I 

(x c + l- k 2 ,T c ) 
i 

x c +l- k 2 ,T c + l) 
i 



It follows from the definitions that 

oo 1 



U U uf^^^^^ 1 '^^^- 

t \cA h, — n j-„— n \ / 



(22) 



(x c ,T c )GAfci=0fc 2 =0 

The weights of introduced sets obey the following relations 

/i(A„(0,0)) = fi(Y(x c ,T c )), 

/x(A,(0,l)) =- /i(W(x c ,T c )), 
/ u(A w (A;,0)) =- n(A v (k + l,0)), k = 0,1,2... , 
/i(A w (A;,l)) =- /i(A„(A; + l,l)), fc = 0,1,2... . 

Using these relations, we can express the probability 

P T (x \x°) = m§ 2 \ u 2 ) 

via the weight of all unallowed configurations 

OO / 

MU 2 )= ^(/^(A,(A;,0))+/i(A w (A;,0))-/i(A,(A;,l))-/i(A w (A;,l))). 

(x c ,T c ) k=0 



(23) 
(24) 
(25) 
(26) 

(27) 
(28) 
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Our next aim is to bring this expression to a determinant form similar to (|15|). Consider 
operator Sj, which acts on the set of the free trajectories of i-th particle and gives the set of 
free trajectories with the origin shifted by 1 in negative direction: 

aiT T (xj\x1) = T T (xj\xi - 1) (29) 

The operator representation allows us to write (}28|) in the compact form: 

MU 2 ) = // {jZ^-Mx2\x° x ) ® T T ( Xl \x 2 )^ , (30) 

where the denominator is defined by its expansion 

1 oo 

T^ = I> < 31 > 

and the summation of operators means the joining of the sets. 

The crucial property of the operator expression in Eq. ()30|) is its factorization with respect 
to indices 1 and 2. Introducing the functions 

oo 

F 1 (x 2 - x°\T) = fi ((1 - aO^TrCxalx?)) = ^ F (x 2 - x° + k\T), (32) 

k=0 

and 

i 

F_i(zi - x°|T) = /i ((1 - a 2 )T T { Xl \x 2 )) = ^(-l) fc F (x! - x° 2 + k\T) (33) 

k=0 

we obtain the two-particle probability Pt(x\x°) = fi(S 2 ) — /i(U 2 ) in the form 

P T (x\x°) = F (xi - x\\T)F {x 2 - x%\T) - F x {x 2 - x\\T)F- X {x x - x° 2 \T) (34) 

or 

p T (x = rfet M, (35) 

where 

M y = (xj - x°|T) z, j = 1,2 (36) 

To generalize the determinant formula (J35j) to the general P-particle case, we should take 
into consideration several additional properties of intersecting trajectories. First, we may 
organize the procedure of exclusion of unallowed elements of Up in an ordered way. Starting 
with the top line of the lattice A, we examine all free trajectories of the set Sp, row by 
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row, until we meet the first collision point where unallowed trajectories are cancelled with 
elements of the auxiliary set A 7r ( 12 ...p). 

If the number of particles P > 2, the elementary squares associated with the collision 
points (x c , T c ), [(x c , T c ), (x c + 1, T c ), (x c , T c + 1), (x c + 1, T c + 1)] for different pairs of trajec- 
tories may occur several times in one horizontal strip of A. If squares filled by interacting 
trajectories are separated one from another by a gap of empty sites, the above arguments 
can be applied to each pair of interacting trajectories independently. The crucial case for 
the Bethe ansatz is a situation, when the elementary squares are nearest neighbors. The 
specific property of the totally ASEP is that, in each pair of interacting trajectories, the 
right trajectory remains free and interacts with the next trajectory independently on its left 
neighbors. Therefore, we can analyse the interaction between particles considering succes- 
sively elementary squares in each row from left to right starting from an arbitrary empty 
square until all unwanted trajectories will be removed. 

After T' steps from the top to bottom, one obtains the set of path configurations q G Sp, 
which are allowed in the first T' rows, and the set of auxiliary configurations yet not involved 
into cancellation procedures. Remembering that all elements of A 7r n 2 ...p) for all tt ^ 1 have 
end points permuted with respect to the original order 1, 2, . . . , P, we conclude that each 
element of A 7r ( 12 ...p) contains at least one collision point. Therefore, all elements of the 
auxiliary set will be cancelled after T' — T steps with unallowed elements of Up. 

The described way of exclusion of unallowed configurations implies a successive construc- 
tion of the auxiliary set A 7r (i 2 ...p). We notice that each two intersecting trajectories in A 2i 
are non-equivalent: one of them belongs to the particle which overtakes another and can 
be called "active". On the contrary, the second particle is "passive". In the case P > 2, 
one trajectory can overtake m others, and we call it "m-active". Similarly, the "m-passive" 
trajectories appear. 

Assume, that the trajectory of a given particle has m active intersections. It means that 
it participates m times in the cancellation procedure and its starting point is shifted m times 
to arbitrary distances in the negative direction. As a result, the auxiliary set associated with 
the free trajectory between points x® and Xj becomes 

7- \ , MxjVi) (37) 
(1 -a,-) m 
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Similarly, for trajectories having m passive intersections we get 

(l-a 4 ) m r T (x J |x°) (38) 

Note that these expressions may be combined into one, if we assume that m > for active 
trajectories and m < for passive and m = for free ones. The weights of sets of m-active 
and m-passive trajectories are given by functions introduced in jl^ : 



F m (xi ~ 4\T) = jt( k+ m m _ 1 *Wi -x° + k\T) 

h—n \ / 



(39) 



fc=0 

for m > 0, and 



F m (x i -x?|T) = ^(-l)*^ ™y Q { Xj -x<l + k\T) (40) 

for m < 0. Activity m of each trajectory is defined uniquely by the permutation tt(12 . . . P), 
so we have for the weight of auxiliary set A 7r (i2...p) 

p 

M(A 7r (i2...p)) = Y\_ F*(i)-i{xir{i) - x°i\T) (41) 

8=1 

This product together with sgn(ir) is a term of expansion of the determinant det M with 
matrix elements for any permutation differing from the identical one 

M itj = Fi-j(xi - x°j\T) i,j = 1, 2, . . . , P (42) 

It follows from Eq. © that the term corresponding to the identical permutation ix = 1 is 

p 

n(S P ) = Y[F (x i -x° i \T) (43) 
i=i 

Collecting the contributions from the set Sp and all auxiliary sets A 7r ( 12 ...p), we obtain the 
the determinant formula 

P T (x\x^) = det M, (44) 



which is valid for al 
been considered in 



P > 1. The TASEP with equal hopping probabilities on the ring has 



We skip the details of this derivation and consider ring geometry 



just for the TASEP with particle-dependent hopping probabilities. 
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III. THE TOTALLY ASYMMETRIC EXCLUSION PROCESS WITH 
PARTICLE-DEPENDENT HOPPING PROBABILITIES 

A method of solution of the TASEP with particle-dependent hopping probabilities follows 
in many details that of the previous section. As above, we start with the case of two particles 
P = 2 on the infinite lattice. The conditional probability (}2*Tj) is expressed as the weight of 
all allowed subsets /i(§ 2 \ U 2 ), which depends on the hopping probabilities of two particles 
t>i, t>2- The set of all possible free trajectories has the weight 

At(S 2 ) = Qx ||S 2 ||, (45) 

where 

Q=f[vr x hi-vi) T - xi+x °. (46) 

i=l 

The weight of unallowed set /x(U 2 ) is determined by the structure of subsets Y(x c , T c ) C U 2 , 
W(a; c ,T c ) C U 2 and auxiliary sets A v (ki, k 2 ) C A 21 and A w (ki,k 2 ) C A 21 defined by Eqs. 
(ED); (J2~2~|) (see also the text between these expressions). The idea of evaluation of /u(U 2 ) is 
to extract the factor Q from the weight of U 2 , reducing the problem of evaluation of /i(U 2 ) 
to consideration of the weights of bonds at the collision point (x c , T c ) only. 

Extracting Q from the weights of free trajectories, we get unweighted skeletons of tra- 
jectories, with fixed starting and ending points. Thus, we may consider a cancellation of 
skeletons of unallowed sets of trajectories instead of the cancellation of weighted ones. This 
is true also for subsets V(x c ,T c ) and A„(A;i, k 2 ) as they consist of free trajectories weighted 
in a standard way: the step of i-th particle has weight Vi, the stay of i-th particle has weight 
(1 — Vi). However, the step [(x c ,T), (x c ,T + 1)] in W(x c ,T c ) and A w (0, 0) having originally 
the weight — V{, becomes — after extraction of Q, which leads to a correction factor 
in the total weights of sets W(x c ,T c ) and A w (0,0). Therefore, the reduced weights of sets 
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remaining after extraction can be written as 



£(§2) = INI (47) 

£(V(x c ,T c )) = ||V(z c ,r c )|| (48) 

/2(A„(0,0)) = 11^(0,0)11 (49) 

A(W(x c ,T c )) = - (7^7) l|W(x c ,T c )|| (50) 

A(A m (0,0)) = - ( 7^-^ 11^(0,0)11, (51) 



1 — 1> i 

where ft (■) is defined by the identity /1 (■) = Q fx (■) . Along with reduced weights of A„(0, 0) 
and A^(0, 0), we define fx (A 1) (fc 1 , k 2 )) and fx (A U) (fc 1 , fc 2 )) for = 0,1,2... and k\ = 0, 1: 



1 - Ui 



/2 (4,(0,1)) = /2(W(x c ,T c )), (52) 
ft (A w (k, 0)) = /i (A„(A; + 1, 0)) , k = 0, 1, 2 . . . , (53) 

;| ( /i (A W (A;, 1)) = fx (A v (k + 1, 1)) , k = 0, 1, 2 . . . . (54) 

Using these relations, we can write the weight of unallowed set in an operator form and 
verify the identity 

MU 2 ) = Q x J \l ^2 r T (x 2 \x°) ®T T { Xl \x° 2 ) ) (55) 

An important property of Eq. (|B3j) is that the reduced weight of two interacting trajectories 
depends on the weight of only one of them, namely on that named "active" in the previous 
section. Therefore, given a permutation of P end points we can evaluate easily the reduced 
weight of each trajectory taking into account its all active and passive intersections with 
other trajectories. Another important property of Eq. (jHHjl is that the operators a\ and d 2 
enter in a factorized form. This property guarantees that the full reduced weight equals to 
a product of the reduced weights of all trajectories. Fig. [T] illustrates the rule of evaluation 
of reduced weights. 

Consider the first trajectory from x\ to £4. It has three active intersections with tra- 
jectories starting at x 2 , x® and x\. After each collision of two particles, only allowed parts 
of their trajectories survive so that the hopping probabilities V\, v 2 , V3, V4 are ascribed to 
different pieces of trajectories as shown in Fig. ^ Therefore the first trajectory has the first 
intersection with weight depending on i> l3 the second one on v 2 and the third on t> 3 . 
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FIG. 1: The set of four trajectories, which are active or passive in different intersections. 

The second trajectory from x® to 23 has one passive and two active intersections so 
its total activity is +1. The passive intersection has weight depending on v±, which is 
compensated by the weight of the first active intersection also depending on v\. The weight 
of the second active intersection depends on t> 2 . Two other trajectories starting at x® and 
x° have two passive intersections each, depending on v±, v<i and V2, v% correspondingly. We 
may conclude that the reduced weight of each trajectory is defined by its number and its 
activity. The weight of z-th m-active trajectory (m > 0) depends on Vi and m — 1 successive 
values Vi + i, . . . ,t>j +m _i. The weight of i-th passive trajectory of activity — m depends on 
t>;_i, v j_ 2 . . . , fj_ m . Then we may write the resulting expressions for the reduced weights in 
the operator form: 

F itm {xj - xlT) = £ m (l ~ 1 T v k a MxjlXi) J , m>0, (56) 

V i,—n \ ^i+k J ] 



>fc=0 
-1 



Fi, m (xj - xl T) = jl ( TT (l - - Vi+k a,) MxjlxV) ) , m < 



(57) 



and 



Fi, {xj - xl T)=fl {T T { Xj \x°)) . (58) 
Analytical expressions for F it7n (xj — x®,T) are 



F i>m (xj-xlT) = f + z*^ll (l--J^±_^ m>0, (59) 



fc=0 



2niz ' 



Fv.b-w-f H) r ^n O-^)^ m<0 



(60) 
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and 




Here the integration is over the circle around origin with vanishing radius. Computing the 
weights for all permutations, one gets the P x P matrix M with elements 

M itj = F hJ _ t ( Xj - xl T) (62) 

and the conditional probability 

P T (x |f°) = /i(§ P \ Up) = Q x /2(§p \ Up), (63) 
where the weight factor Q generalizes Eqs. (fTBj) and (J4*6*|) : 

Q = f[v?~ X hl-v % ) T -^. (64) 
i=i 

Finally, the explicit expression for Pt(x\x°) generalizing Eq. PHI is 

p 

P T (x\x°) = Y[v. l - X *{l - Vi ) T - Xi+x ° detM. (65) 
-i=i 

In the continuous time limit (T — > oo, Vi — > 0, < w^T < oo, u^T ~i G M + , i = 1, 2 ... P), 
we come to the result obtained by Rakos and Schiitz [l|. 

To extend the derivation to the ring of length L, we map the trajectories wrapping the 
cylinder A on the infinite integer strip Q = [— oo, oo] x [0, T] introducing extended coordinates 
X = x + nL, n integer, for equivalent points. The initial coordinates of particles are x° = 
(xi, x®, • • • , x° P ), < X® < < . . . < x*p < L — 1. The positions of particles on the ring at 
time T are given by ordered coordinates x = (x\, X2, ■ ■ ■ , xp), < X\ < X2 < ■ ■ ■ < xp < L—l 
and cyclic permutation of indexes of particles a = (a%, a 2 , . . . , ap). The ring coordinate of 
i-th particle {i — 1, 2, . . . , P) at time T is x ai . Let z-th particle intersects the bond [L— 1, 0] 
n* . times during time T. Then the number of steps it advanced is 

N, = X ai - xl (66) 

where X ai = x ai + n* a .L is the extended coordinate of i-th particle on the strip. Since the 
particles cannot overtake one another, we have for any cyclic permutation a: 



X ai < X a2 < . . . < X ap < X ai + L. 



(67) 
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Clearly, these inequalities impose restrictions on rotation numbers n* = (n*,^, . . . ,n P ). 
Particularly, the set of independent variables describing unambiguously the state of system 
are the ordered ring coordinates x and the total current through bond [L — 1, 0] 

p 

J = $>*. (68) 

i=l 

Our aim is to find the conditional probability to find the particles in the ring positions for 
a fixed J at time T if the particles start from initial positions x°. 

Consider a set of configurations of P free trajectories on the strip Q together with their 
copies placed periodically with period L in horizontal direction. We obtain the set of con- 
figurations 

§ P = T T (x ai + n* ai L\xl) ® T T {x a2 + n* a2 L\x° 2 ) <g> . . . <g> T T {x ap + n* ap L\x° P ). (69) 

repeated on the strip infinitely many times. As above, a configuration is allowed if there are 
no collision points among trajectories. For an allowed configuration, there is at least one 
possibility to draw it without collisions. To eliminate the unallowed set Up of free config- 
urations, we introduce a new auxiliary set. In the case of periodical boundary conditions, 
this set includes not only the permutations of numbers of particles tt, but also all possible 
numbers of rotations n = (n 1; n 2 , . . . , np): 

K(i,2,-,P),n = T T (x 7r{1) +n 7r{1) L\x 1 )®T T (x 7r{2) +n 7T{ 2)L\x 2)'S)- ■ .®T T {x^ P) +n^ P )L\x P ) (70) 

The numbers of rotations are arbitrary but yield the condition 

p 

J = 5>*> ( 71 ) 

i=i 

because the permutation of end points in any pair of trajectories does not change the sum 
of rotation numbers. The weights of configurations Sp and ^(1,2,. ..,p), a a s well as their 
reduced weights can be defined in a way quite similar to that for the infinite lattice (Fig. 
12} . The difference between the infinite lattice and the ring is a possibility for given particle 
to overtake another one several times. For instance, the trajectory x\ — > £2 in Fig. [2] has 
two active intersections with both trajectories from x® to x± and from x® to 23. So its 
total activity equals m — 4. The trajectory x° — > £3 has three passive intersections and its 
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FIG. 2: The configuration of three trajectories on the ring. 

activity is m = —3. The trajectory has two passive and one active intersection, so 

its activity is m = — 1. In general, the activity of trajectory xf — > Xj is given by 

p 

Sij =j~i + Prij - ^ n k- ( 72 ) 

k=l 

The weights of different pieces of the trajectory starting from x° vary from to the next 
values Vi+i or Vi-i for each active or passive intersection. If % = P, the next value Vi+i = v± 
and if i — 1, the value t>j_i = tp. Generally, the indexes of hopping probabilities coincide 
by modulo P: 

Wi +fc P =Vi, k = ±1, ±2, . . . (73) 

The sign of a configuration is defined by the sign of permutation of end points. Using these 
rules, we derive the reduced weight of allowed trajectories in the form 

/2(Sp \ Up) = /2(S P ) + ^ ^(-l) 1 ^ 1 ' sgn(n) A(A 7r(12 ...p ))fJ ), (74) 

where the first summation is over all integers {rij £ Z; i = 1, 2, . . . , P} that satisfy the rela- 
tion (j71jl . The second sum is over all permutations 7r of particles except initial permutation 
<x This expression can be written as a sum of determinants. Introducing the matrix M 

M id = F itSi . ( Xj + njL - xl T) , (75) 

we obtain 

fi(Sp \ Up) = i- 1 )^ det M - ( 76 ) 
ft 
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The transition probability for fixed J is defined by identity 

P T (x, J|x°) = n(S P \ Up) = Q x ft(Sp \ Up), (77) 

with the factor 

Q = f[v?'{l-Vi) T - Nt - (78) 
i=i 

Finally, we find 

P T In -n I 

P r (f, J|f°) = ]Jv^(l - Vi) T - Ni Y^(-l) i<j ' ™ J det M. (79) 

Being calculated explicitly, Eq. [7U] is a polynomial in i>x, . . . , Vp. For small L,P,T, it can 
be found directly from consideration of allowed configurations of trajectories and compared 
with results of evaluation of Eq. EH For example, in the case L = 5, P = 3, T = 4, 
x° = (0, 1, 2), x = (0, 1, 2) and various J, we get 

P T (f,J = 0|f°) = (l-t;3) 4 , (80) 



Fr(x,J=l\&)= v lV2 v 3 3 (1 - v 2 f {vfv 2 + vfv 3 - 2v\ - 7v\v 2 + (81) 
2vfv 2 v 3 - 7v\v% + \2v\ + 18viV 2 - SvxV 2 v 3 + 
18v!V 3 - 28^i - 22v 2 + 12v 2 v 3 - 22t> 3 + 32) , 



P T (x, J = 2|f°) = 6(1 -vxfvjvfrl (82) 

For J > 2, the probability Pt{x, J\x°) vanishes. 

Summation over J gives the transition probability from an arbitrary initial ring coordi- 
nates x~o to final coordinates x: 

oo 

p T (f = Pt{x, A^)- (83) 

J=0 

In the specific case, when the hopping probabilities Vi are not particle-dependent, we can 
write Q in the form 

p p 

V (xi-x°+n*L) T- V (xi-x°+n*L) 

Q = v i=1 (1 - v) 1=1 . (84) 
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The conservation condition 

p p 

i=l i=l 

makes it possible to change all variables n* by m and evaluate the sum over J. After the 

nn 

summation, we recover the known result for equal hopping probabilities |15L 1161]: 

+oo +oo V I — ^ 

P T (f|f°)= ••• E ( _1 ) 1<J ^ ^ Y[v x *- x > +n * L (l -v) T -^ +x °^ L detM. (86) 

n\=— oo np=-oo i=l 

For considered values of parameters, it gives 

p T (x = (1 - w) 2 (6w 10 + 4t> 9 - 24v 8 + 60v 7 - 72v 6 + 32v 5 + v 2 - 2v + l) . (87) 



IV. DISCUSSION 



We have calculated exactly the non- stationary conditional probability for the TASEP 



with particle-dependent 
been studied by Evans 



lopping probabilities on a ring. The steady state of this model has 
, who found that the mean velocity v s t in the thermodynamic limit 



(T — > oo, L — ► oo, P — > oo, P/L = p) is determined from equation 

1 1 — Vi , , 

1 -p = ^tE 88 

£ - v st 

i=i 

Obviously, v st does not depend on the permutation of velocities v = (vi, t>2, . . . , vp). In the 
low-density phase the fast particles are sticking to slower ones and form a number of jams just 
behind slow particles. The slow particles are mainly free, i.e. almost all of them have empty 
neighboring sites in direction of motion. In the thermodynamic limit, the jam corresponding 
the slowest particle dominates over other clusters. This behaviour has a strong analogy with 
Bose condensation and the steady-state velocity of a particle is equivalent to the fugacity of 
an ideal Bose gas 0. For high densities, the slowest particle is not free, the jam dissolves 
and a congested phase appears. 

The non-stationary state of the TASEP considered in this paper differs from the station- 
ary case in several respects. First, the non- stationary mean velocity depends not only on 
absolute values of particle velocities, but on their order too. Second, a time dependence of 
the total number of particles in the jam is an important characteristic of the traffic 19. 2p|. 



Generally, the exact non-stationary solution of the TASEP on a ring may describe how the 
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phase transition arises in thermodynamic limit. A next possible step in this study would be 
the TASEP with parallel dynamics, which answer the question whether distinct updating 
schemes produce different collective behaviour. 
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